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Abstract 

We give a simple and entirely elementary proof of Gasper's theorem on the Markov sequence 
problem for Jacobi polynomials. It is based on the spectral analysis of an operator that arises 
in the study of a probabilistic model of colliding molecules introduced by Marc Kac, and the 
methods developed here yield new estimates relevant to the collision model. 

Mathematics Subject Classification Numbers: 31B10, 33C45, 37A40 

1 Introduction 



O ' 1.1 The Markov Sequence Problem and the Theorems of Bochner and Gasper 

oo ; 

' Let (X, be a probability space. A Markov operator T on L'^ifJ.) is a linear operator that 

preserves positivity; i.e., / > =^ Tf > 0, and preserves the constants; i.e., Tl = 1. If T is self 
^ ■ adjoint, it follows by duality and the Riesz-Thorin interpolation theorem that T is a contraction 

on LP(^) for all 1 < p < oo,) Consequently, the spectrum of T lies in the interval [—1, 1]. 

The next definitions, which are less standard, are taken from [3] and [Ij: A unit orthonormal 
basis for is an orthonormal basis {fn}n>o such that /o = 1. Though we discuss a broader 

class of examples in Section 5 and in the Appendix, In the main examples here, X = M, or some 
subset of M, and {/n}n>o is the sequence of orthonormal polynomials for ^u. In any case, we shall 
always suppose that X is a locally compact Hausdorfi^ space, and that is a Borel measure. 

Given a unit orthonormal basis {fn}n>o, the set of Markov sequences Ai for this basis is the 
set of all sequences {A„}n>o such that there exists a self adjoint Markov operator K with 

Kfj = Xjfj for ah j>0 . 
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Notice that necessarily Aq = 1 and A„ G [—1, 1] for all n. Also, since a convex combination of self 
adjoint Markov operators is self adjoint and Markov, Ai is convex, so that Ai may be described 
by specifying its extreme points. 

The Markov sequence problem is to determine, for a given unit orthonormal basis, the set Ai. 
Naturally, it is sufficient to find the extreme points. 

The Markov sequence problem seems to have been first considered by Bochner [6], and the first 
result, for ultraspherical polynomials, is his as well. 

We recall that the for each 7 > —1/2, the ultraspherical polynomials {pn^}n>o are the orthonor- 
mal polynomials, for the measure ^^"''^ 

d/x(^)(t) = c^(l-t2)T-i/2dt where = ^ ^ (1-1) 

^ ^ ^ ^ ^ V^r(7 + i/2) ^ ^ 



is the normalization constant that makes probability measure. The normalization as unit 

vectors in L'^{fi^"'^) is just one useful and frequently encountered normalization. Another that will 

(7) 

be useful here is generally denoted with an upper-case P: The ultraspherical polynomials Pn are 
normalized so that pjC\l) = 1 i.e. 

P^^'H-) = ^- (1-2) 

Pn [l) 

Throughout the paper, an upper-case P denotes this normalization, while a lower case p denotes 
the L?'{ix^'^^) normalization. 

The ultraspherical polynomials are special cases in the wider family of Jacobi polynomials: 
Recall that the Jacobi polynomials p"n'^^ form an orthonormal basis for I? ([—1, 1], dfi"'^^ where 

/,"'/3(dx) = c«,^(l - + xfdx , (1.3) 

where Ca^p makes probability measure. In particular, the ultraspherical polynomials arise 

for the special case 7 = a — 1/2 = ^ — 1/2; that is 

P^r:'\t) = p^y-y^'^-y^Ht) . (1.4) 

Theorem 2 of [6] may be phrased as follows: 

1.1 THEOREM (Bochner). For any 7 > 0, the sequence {A„}„>o is a Markov sequence for 
{Pn^}n>o if o-nd only if there is a probability measure u on [—1,1] such that 

' P'^^^^huit) . (1.5) 



ipi7\i) 

For each such Markov sequence {Xn}n>o, the measure v is unique. In other words, for each t, 
{Pn\t) /pn\^)}n>o is a Markov sequence for {pn^}n>o, CLnd these are the extreme points of the set 
A4 of all such Markov sequences. 

Since the ultraspherical polynomials are Jacobi polynomials with a = /3, it is natural to ask 
whether one can one extend Bochner's result to a wider class of Jacobi polynomials with a ^ (3. 
This question was answered by Gasper pTj [T2]: 
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1.2 THEOREM (Gasper). For a > 13 with 13 > -1/2 or a > (3 with 13 = -1/2 the sequence 
{A„}n>o is a Markov sequence for {pt^^} n>o, if and only if there is a probability measure v on 
[—1,1] such that 

For each such Markov sequence {A„}„>o, the measure v is unique. In other words, for each t, 
{pr^'^\t) /p^n'^\l)}n>Q is a Markov sequence for {pn^}n>o, a,nd these are the extreme points of the 
set M of all such Markov sequences. 

While Bochner's proof of Theorem ll.ll is not terribly complicated, Gasper's proof of Theorem 1 1.2 1 
is far from elementary. Even though it has been simplified by the work of others, particularly 
Koornwinder, it remains a tour de force: Koornwinder's proof still uses many deep results on 
special functions. 

In this paper we shall give entirely elementary and self-contained proofs of these theorems. 
Moreover, these proofs will allow us to obtain bounds on the sizes of the extremal eigenvalues. 
Before stating the new results more precisely, we recall the proof of Bochner's Theorem, as this 
will clarify the matter of what was already well understood, and what was in need of clarification. 

1.2 Product formulas and the Markov sequence problem 

In this subsection we explain that the Markov sequence problem is easily solved for unit orthonormal 
sequences that satisfy a product formula, as defined below. Indeed, Bochner's original proof of his 
theorem went by this route, and was facilitated by the fact that the product formula he required 
had already been established long ago by Gegenbauer. Gasper, on the other hand, had more work 
to do since before his work, no general product formula for Jacobi polynomials was known. 

1.3 DEFINITION (Product Formula). A unit orthonormal sequence {/n}n>o satisifes a product 
formula in case there exists function {x,y) i— )■ d^x^yiz) from X x X to the space of probability 
measures on X, and also some xq £ X such that for each n > 0, 



Fn{x)Fn{y) = / Fn{z)dfXx,y{z) , (1.7) 

Jx 



where Fn{x) — "^"^ 



fn{xo) 

For example, in the case of the ultraspherical polynomials {p^n^}n>o, take xq = 1, so that F„ 
becomes PjC^ Then one has Gegenbauer's identity |13j . which dates back to 1875: 

1.4 THEOREM (Gegenbauer's Identity). For all 7 > 1/2, and all n>0, and all a G (-1, 1), 
P;iia)P2it) = £^ P2 {at + sVl^Vl^) d/i(^-i/2)(^^ _ (18) 

To see this as a concrete instance of the abstract product formula (jl.7p . let 5„ denote the Dirac 
mass at u G [— 1, 1], and define 

dMa,(^) = 6^,^,^^{z)dt,^^'-"^){s) . 
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Then ([HD becomes P2{a)P2{t) = j P2{z)diia,t{z) , as in ([TT]). 

The following theorem relates the Markov sequence problem to the problem of establishing a 
product formula. The theorem summarizes ideas that can be found, reading between the lines, in 
Bochner's paper [6] for the ultraspherical polynomials, and much more explicitly, and in general, 
in the paper [3] of Bakry and Huet. 

1.5 THEOREM (Markov Sequences and Product Formulae). Let X be a closed interval 
in M, and let ^ be a regular Borel probability measure whose support is X. Let {/n}n>o o unit 
orthonormal basis for L'^{^) consisting of real valued functions. Let xq be any fixed point in X. 
Then the following are equivalent: 

fJx) 

(1) For each x £ X, {A„(x)}„>o is a Markov sequence for {/n}n>o where Xn(x) := -r^, — r. 

fn{XQ) 

(2) For each {Xn}n>o G there exists a a Borel probability measure u so that 

An = / ^Au{x) ■ (1.9) 

Jx fn{Xo) 

(3) With Fn{x) := fn{x)/ fnixo), the {Fn}n>o satisfy the product formula d^.Tp for some family 
d^x,y{z) of probability measures on X. 

Finally, if any (and hence all) of these conditions are satisfied, and if {fn}n>o is a sequence 
of bounded continuous functions whose finite linear combinations are dense in Cb{X), then the 
probability measure v in \1.9^ is unique, so that M. is a simplex and the {fnix) / fn{xo)}n>o o'^e its 
extreme points. 

The equivalence of (1) and (2), as well as the statement concerning uniqueness of the measure z^, 
is due to Bakry and Huet, [3], together with many other results on the Markov sequence problem. 
The equivalence of (1) and (3) is implicit in Bocher's paper |6|, though his argument is different 
from what follows below, and in particular, he makes no use of self-adjointness of certain operators 
associated to product formulae - a crucial feature of our approach. Thus, while we make no claim 
of originality for the results in Theorem II. 5( we provide a complete proof for completeness and 
clarity. 

Proof: We first show that (1) implies (2): If {fn{z) / fnixo)}n>o is in Ai, then by the spectral 



theorem, Kz{x,y) = A„(z)/„(x)/„(y) is the kernel of a Markov operator on L'^{p) with 

n=0 

Kzfn{x) = Xn{z)fn{x). If is a Borel probability measure then K := Kzdv{z) is a Markov 

Jx 

operator with eigenvalues 

/ \n{z)Mz) = [ jj^duiz) . (1.10) 
Jx Jx fn{xo) 

We next show that (2) implies (1): Suppose that {A„}„>o is a Markov sequence. Again, by the 

oo 

spectral theorem K{x,y) := Xnfn{x)fn{y) is the kernel of a Markov operator K on L'^{fi) with 

n=0 

Kfn{x) = \nfn{x). Since K is Markov, K{x,y) > for all x and y, and for each y, K(x,y)dfj,{x) 
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is a probability measure. Taking y = xq, define the probability measure dv = K{x, xo)dn. Then, 
for each k, 



/ fk{x)diy{x) = fk{ 

JX JX 

oo 

= E 



^ A„/„(2;)/„(xo)d/i 



.n=0 



n=0 



fk{x)fn{x)dll 



Kfn{xo) = Xkfk{xo) , 



(1.11) 



which proves (|1.9p . 

We next show that (1) implies (3): If {/n(^)//n(a^o)}n>o is in Ai, then for each z in X, 



kz{x,y) = J2 



n=0 



fn{z)fn{x)fn{y) 
fn{xo) 



is the eigenfunction expansion of the kernel of a self adjoint Markov operator. Evidently, (x, y, z) i— )■ 
kz{x,y) is invariant under any permutation of x, y, and z. Thus, for each m. 



J^fm{z)h{x,y)dfi{z) = J^k^{y,z)fm{z)dfi{z) = J^J^^f^^y^ ■ 



(1.12) 



fm{xo)' 

To recognize this as a product formula, for each x, y, define a probability measure fix.y by 

diJ,x,y{z) = kz{x,y)dti{z) . (1.13) 

Then with Fn{x) := fn{x)/ fn{xo), (frT2]l becomes ([LTD . 

We next show that (3) implies (1): Assuming (3), fix y G and define an operator Ky on Cb{X) 

by 

Ky(p{x) = / (p{z)dfix,y{z) . 

JX 

Since K is a Markov operator, it has a bounded extension to L'^{fi). The product formula says 
that for each m, fm is an eigenfunction of Ky with eigenvalue Fm{y)- Any bounded operator with 
a complete orthonormal basis of eigenfunctions, each of whose eigenvalues is real, is necessarily self 
adjoint. Thus Ky is a self adjoint Markov operator, and hence {-Fn}n>o belongs to M. 

Finally, It remains to show that the measure v is uniquely determined. For this, let / be any 
continuous function on X. Let e > be given, and let g{x) = 'Ylin=o'^nfn{x) be a finite linear 
combination of the {/n(a^)//n(a^o)}n>o such that |/(x) — g{x)\ < e for all x. 

Let and d be two Borel probabilty measures such that (jl.9p holds for some {Xn}n>o G Ai. 
Then 

» N N N 

/ g{x)du{x) = ^an /„(x)di^(x) = ^ a„A„/„(xo) = ^ a„ / fn{x)du{x) = g{x)dV{x) . 

•'^ n=0 •'^ n=0 n=0 •' ^ •' ^ 



Therefore, | f{x)dv{x) — J-^ f{x)d'D{x)\ < 2e. Since e is arbitrary, j-^ f{x)du{x) = J-^ f{x)d'D{x) 
for all / G Cb{X). This of course means that = P. ■ 
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1.6 REMARK. Notice that by the Weierstrass Approximation Theorem, the conditions in the 
final part of Theorem \1.5\ are automatically satisfied in any application to orthogonal polynomials 
on a compact interval. 

It is interesting to note that once one has the product formula (|1.7p . one can use it to define a 
convolution: For any two finite, positive measures A, on X, define the convolution A * of A and 
u by 

A*z.:= J dfi^,yiz)dXix)diy{y) ; (1.14) 

this too is a finite positive measure. 

Note that the "Fourier" coefficients of A and z^, given by J F„dA and J F„dz^ satisfiy 



J FndX J Fndu = J Fnd{\ 



so that the usual relation between Fourier coefficients and convolutions holds. We now return to 
the matter of proving Bochner's Theorem. 

1.3 Gegenbauer's Identity and a proof of Bochner's Theorem 

Proof of Bochner's Theorem: Since we have a product formula for the ultraspherical polyno- 
mials, namely Gegenbauer's identity (II. Sp . condition (3) of Theorem [L5] is satisfied, and Bochner's 
Theorem follows immediately. ■ 

This proof is simple, but hardly complete: one needs the product formula. In the case of the 
ultraspherical polynomials, this was ready at hand since 1875. For the general case of the Jacobi 
polynomials, no product formula was available when Gasper began his work. His strategy was to 
show that for {/n}n>o being a sequence of Jacobi polynomials, one has the positivity result 

ff^mpfiM>o, (1.15) 

holding pointwise almost everywhere. Then, this sum defines a positive kernel, which can be used, 
as in the proof of Theorem [L5l to prove a product formula. However, this direct proof of pointwise 
positivity is far from simple. 

In this paper we present a truly simple approach to the product formula for Jacobi polynomials. 
First, however, we shall illustrate this approach by providing a simple, self-contained proof of 
Gegenbauer's identity (|1.8p . 

For our purposes, it is most helpful to consider (|1.8p as an eigenvalue identity. 

1.7 DEFINITION (The Correlation Operators). For each 7 > 0, and a G (-1,1), define an 
operator Ka on L'^{ij,^"'^) by 

Kaf{t) = f (at + sVl-aVl-t^) d/x(^-i/2) _ 
We refer to the Ka as the correlation operators for reasons that will be explained in Section [2j 
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With this definition and (|1.2p . ()1.8|) can be written as 

KaP^;;'\x) = ?^p(;y\x). (1.17) 

Pn (1) 

Thus {p^7^(a)/p^7^(l)}n>o is the eigenvalue sequence of Ka, and since the eigenvalues are real and 
the eigenfunctions are orthogonal, it follows that Ka is self adjoint, as noted in general in the proof 
of Theorem 11.51 

Proof of Gegenbauer's identity: The starting point is a direct proof that Ka is self adjoint. 
Prom (I1.16P we find, 

{Kaf, = c^c^_i/2 g{t)f [at + sVl-a^VT^) (1 - s^r^'dsil - t^y-'/^dt . 

With the change of variables u = at + sVl — a^\/l — t^, the integral over s becomes 

at-VT^VT^ (l-a^)T 
so that 

.^ , ,((l-a2)-(n2 + t2_2atu))X"^ , , 

{Kaf,g)L^^^(-,)) =C^C^-l/2 y J ^9{t)f{u) (1 _ a2)7-l/2 ' 

where (•)+ denotes the positive part. Thus, Ka is self adjoint on L'^{fi^'^^). 

For the rest, we reurn to the original defining formula (11.160 First, by the symmetry properties 
of Ka maps polynomials of degree n to polynomials of degree n. It follows that the 

spectrum is discrete and the eigenfunctions are polynomials that are orthogonal with respect to 
the measure IjL^'^\ and hence are the pn^ . Let A„ be the eigenvalue corresponding to pn^; i.e., 
^nPri\t) = KaPn\t)- Taking the limit t — )• 1 on both sides, using the Dominated Convergence 
Theorem yields, 

A„P^7Hl)=P^7Ha) (1.19) 

which immediately gives Gegenbauer's identity in the form (|1.17p . ■ 
In what follows, we shall make repeated use of the mechanism illustrated in our proof of Gegen- 
bauer's identity, and the next theorem paves the way for its broader application: 

1.8 THEOREM (Evaluation Formula). Let X be a closed interval in M, and let be a regular 
Borel probability measure whose support is X. Let {/n}n>o tl^^ "anit orthonormal basis for L'^{fi) 
consisting of the normalized orthogonal polynomial for fi. Suppose that for each z £ X , K^ is an 
operator on L'^{fi) with the following properties: 

(1) Kz is self adjoint on L'^{p). 

(2) Lf f is a polynomial of degree no greater than n, then so is K^f ■ 

(3) There exists an xq £ X such that for any continuous function f , and any z £ X , 

hm KJ{x) = f{z) . (1.20) 

fn{z) 

Then for each n, fnixo) ^ 0, and for each x Kzfn{x) = -j—, — zfnix), so that if K^ is a Markov 
operator, then {fn{z)/ fn{xo)}n>o is a Markov sequence for {fn}n>o- 
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Proof: Properties (1) and (2) immediately imply that each is diagonalized by polynomials that 
are orthogonal in so that the eigenfunctions of Kz are the /„. To determine the eigenvalues, 

start from the definition of the nth eigenvalue A„, Kzfn{x) = ^nfni-c), and take the limit x — )• xq. 
By (3) we obtain fn{z) = XnUixo), which tells us fn{xo) 7^ and An = fniz)/fnixo) . ■ 
To summarize, through an analysis of the operators Ka, based on the three properties high- 
lighted in the previous theorem, we obtain a self-contained proof of Gegenbauer's identity, and 
hence Bochner's Theorem. Is there a analogous family of operators that gives Gasper's theorem? 

1.4 Gasper's Theorem 

For /i G C([-l, 1]) define Ka,o by, 

{Kafih){t)= I I h {1 + t) -l + h^ {I- t)r'^ + 2ahr{l-ff/'^ COS e dmo,,p{r,e) (1.21) 
Jo Jo ^ ' 

where 

0) ^ 2T{a + l) _ 2sa-/3-l 2/3+1 . 2/3 ,n Q 22) 

am^,(s[r,U) ^r(/3 + l/2)r(a-/3)^^ > r sm aardti [i.ZZ) 
is a probability measure. We now have, 

1.9 LEMMA. For all a G (-1,1) and a > /3 > -1/2, the operator Ka,o on C([-l,l]) as defined 
in has the following properties: 

(1) Kafl is self adjoint on L'^{fi'^'^). 

(2) The space of polynomials of any fixed degree is invariant under Kafi- 

(3) For any continuous function h, limit's o^(0 — h{2a^ — 1). 

Proof: Given the explicit formula ()1.2ip . the proof of (2) follows from the form of which 
shows that only even powers of cosO are nonzero when integrating over 0. Part 3 follows from the 
evaluation property (|1.20p and the Dominated Convergence Theorem. It is only (1) that requires 
more work. 

We now use a sequence of variable changes due to Koornwinder [16], but for a different purpose. 
We shall contrast our use of it with Koornwinder's in the final section of the paper, but for now, 
suffice it to say that Koornwinder was not concerned with self-adjointness, which is the issue before 
us. 

Consider hi and /i2 in C([— 1,1]) Then by p.2ip and the change of variables t = 2s^ — 1, 
{hi, Kaflh2) l2(^^^ ^-^ is a constant multiple of 

fl rl fn 

q{hi,h2):= / / /ii(2s2 - l)h2{{2a^s'^ - 1) + 26^(1 - s^y + 4a6rs\/l - cos 9) 
Jo Jo Jo 

X sin2/^0(l -r2)°~/3-V2/5+i(l -s2)°s2/^+Msdrd0. 

We must show that q{hi, /12) = q{h2, hi). 

The first step is to replace (1 — r^)"~^~^ by (1 — r^)"^^ and extend the domain of integration 
in r to (0, 00). The point is that we may then regard the integration over r and 6 as an integration 
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over the upper half plane in M^. Changing to Cartesian coordinates x and y yields 

q{hi,h2) = C hi{2s^ - 1)(1 - s'^Ts'^^^^ 
Jo 

(/•OO j-QO \ 
j J h2{2{cx + as)^ + 2c^y^ -l){l-x^ -y%~'^~^y'^^dydx]ds 

where c = 6(1 — s^)-*^/^. 

The second step is to translate and scale, making the change of variables y' = cy and x' = cx+as. 
Then since (1 — s^)°i/^'^dxdy = 6^^"c^("^''^^)dx'dy', this yields 

q{hi,h2) = /\i (2^2 - 1)^2/^+16-2" 







/OO poo \ 
/ h2{2{x'^ + y'^) - l)(c2 - {x' - asf - y'')l-^-'y'"'dy' dx' ds . 
OO Jo J 

Finally, the third step is to change back to polar coordinates; i.e., make the change of variables 
{x',y') — )• (/0, 0). This yields, making crucial use of + 6^ = 1, 



1 

2 _ 1^c2/3+U-2q 



q{hi,h2)= / hi{2s'-l)s' 

Jo 

(roo rn \ 
J J /i2(2p2-l)(62_52_^2^2apscos(/.)+"'^"V2/'+M(/.dpjds 

= 6-2a f' /■\,(2s2-l)/i2(2p2_l) 

Jo Jo Jo 

X (62 -s^-p^ + 2aps cos cP)^-''-' p'"-^' S^/S+l ^-^2/3 ^^^dpds , 

which finally renders the symmetry manifest. ■ 

Proof of Theorem 11.21 The case a = f3 > —1/2 are contained in Bochner's theorem. For 
a > P > —1/2, Lemma 11.91 implies that the family of operators Kafi, a £ [— Ijl] satisfies the 
conditions of Theorem 11.81 for the unit orthonormal basis {plt'^^}n>o- Then from the conclusion 
of Theorem II. 8^ we may apply Theorem 11.51 to obtain Gasper's Theorem in this case. The case 
a >/?,/? = -1/2 follows since p£'"^(x) = pt'~^^^\2x'^ - 1) ■ 
The remaining mystery at this point is where the operators Ka and Kafi came from. In fact, 
the operator Ka arose naturally in the work [8j on the Kac model [H]. The Kac model is a model 
from mathematical physics for the trend to equilibrium in a gas of N molecules interacting through 
binary collisions. An analysis made in [8] of how the rate of equlibriation depends on for one 
dimensional velocities reduces this issue to the determination of the eigenvalues of the operators 
Ka (jl.lGp . which measure correlations between the different particle's velocities. In the analysis of 
the Kac model for three dimensional velocities the following operator naturally arises, 

ICafiv) = j^f{av+ x/l-a2Vl-b|22/) dv^,N-i{v) , (1-23) 

where 

I Qm{N-l)~l\ 

dl^n^Mv) = ^^SmN-ll (1 - IH')^'"^^-^^-'^/'d. , (1.24) 
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and B is the unit ball in MJ", \S'^^^\ is the surface area of the unit sphere in M'^. The number a 
ranges from —1 to 1, > 2 and m > 1. 

Note the similarity of JCa to the operator Ka defined in (jl.l6p . Of course fCa acts on functions 
on the unit ball B, however there is a direct connection to operators that act on functions of [—1, 1], 
such as Kafl- This follows from the fact that /C^ commutes with rotations and therefore preserves 
the class of radial functions. For h £ C{[—1,1]), define Kafih by 



{Ka,oh){2\v\ 



1) 



where 



fiv) := h{2\v\^-l). 



(1.25) 



This operator is well defined since ICa preserves the class of radial functions. A calculation, which 
we shall make in Section 12.31 shows that for 



a 



(m(iV-2) -2)/2 and /3 = (m-2)/2, (1.26) 

the operators -fCa.o defined in (|1.2ip . and (|1.25p are the same. Thus, at least for the half integral 
values of a and /3 in (jl.26p . the apparently more complicated operator defined in (|1.2ip does 
indeed come from an operator bearing a striking resemblance to the one in Gegenbauer's identity. 
Moreover, while m and N are integers in (ll.23p . once the radial part has been rewritten in the form 
()1.2ip . there is no obstacle to letting a and /3 vary continuously. 

The remarkable thing about this construction of Ka^o is that it only uses one invariant subspace 
of the operators ICa to recover the known results of Gasper and Koornwinder. However, there are 
other invariant suspaces with a direct connection to Jacobi polynomials. In fact, we shall see that 
for each integer i > there is a family of operators Ka/ to which we may apply Theorem 11.81 . 
For example, the following result about Jacobi polynomials is the analog of Gegenbauer's product 
formula (jl.Sp : 

1.10 THEOREM. For all a > f3 > —1/2, and all non negative integers £, 



iPn 



it) 



11 



JO 



E 

3=0 



0^(1 + t) + 6^(1 _ t)r2 + 2a6\/l - f^r cos 9 

i/2 



\hr) 



l-t 
T+t 



Pf^(cos0) 



<^rna,p{r,e) , 



(1.27) 



where b = \/l — as before. 

Note that Gasper's formula appears for the case 
have 



0. As a consequence of this theorem we 



1.11 THEOREM. Consider any a > f3 > -1/2 and any integer I > Let P^^^ is the ultras 



pherical polynomial with the normalization 



(1) = 1. Then for all t G [-1, 1], 



pt'^'\i) 



Jo 



t)r'^ 



+ i 



Vi-t2 



r cos ( 




(ir)^'pf^(cos0) 



dma,l3{r,9) . 
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where drua^p is given by il.22\) . 

The case £ = is a well-known formula of Koornwinder [T5j. It was pointed out to us by an 
anonymous referee of a previous version of this work, that for £ > 0, the formula of Theorem 11.101 
is equivalent (though not in such an obvious way) to a product formula due to Koornwinder and 
Schwartz for orthogonal polynomials on the so-called parabolic biangle; their equivalent formula is 
(3.13) in [H]. 

We state and indicate a direct proof, using our methods, of the product formula for the parabolic 
biangle in the Appendix [XI In this Appendix we also discuss the product formula on the triangle. 

The rest of the paper is organized as follows: In Section [2] we give a geometric and probabilistic 
interpretation of the operators Ka and /C^. This shall explain our reasons for referring to them 
as "correlation operators". It also yields a simple proof of their self adjointness, at least for the 
"geometric" values of 7, a and /3. We then use the rotational invariance of the operators ICa to 
determine a sequence of invariant subspaces for them, indexed by the non-negative integer i, and 
we study the spectrum of the restrictions Ka/ of ICa to these invariant subspaces. Though for i > 0, 
Ka/ is not Markov, Theorem 11.81 is applicable nonetheless: The eigenvalues are again expressible 
in terms of ratios of Jacobi polynomials, and in this enable us to easily prove Theorems 11.101 and 
11.111 for the "geometric" values of 7, a and /3. 

Then, in Section [3l we show how 7, a and (3 may be allowed to vary continuously, and thus 
prove Theorems 11.101 and 11.111 in full generality. 

In Section [5] we use the Laplace formula for ultraspherical polynomials and Theorem (jl.lip to 
obtain sharp bounds on ratios of Jacobi polynomials. That is, we obtain sharp bounds on the 
eigenvalues of the extremal Markov operators, and these bounds give sharp information on the 
operator trace classes to which the extremal Markov operators belong. This information is then 
used to discuss the pointwise convergence properties of the eigenfunction expansions for the kernels 
assoicated with the operators in Bochner's and Gasper's theorem. 

In Section [5l we discuss the history of Bochner's and Gasper's results and finally in the appendix 
we state and outline a proof of the parabolic biangle and triangle polynomial product formula of 
|18j along the lines outlined in this section. 

Acknowledgment: The authors would like to thank Dominique Bakry for illuminating discus- 
sions about his papers with Huet and Mazet, and for suggesting that ideas arising in our work on 
the spectral gap for the Kac model with three dimensional velocities might lead to a self-contained 
proof of the Markov sequence problem for Jacobi polynomials. We also thank an anonymous ref- 
eree of a previous version of our paper for pointing out the connection of Theorem 11.101 and the 
parabolic biangle polynomial product formula of Koornwinder and Schwartz |18j . 

2 The geometric cases 

The proof of the product formula for Jacobi polynomials with a = {m{N — 2) — 2)/2 and /3 = 
(m — 2)/2, where m and N are positive integers as in (jl.26p . is particularly simple because of a 
geometric picture for the correlation operator ICa in these cases. In this section, we shall present 
a complete proof for these cases, which we call the geometric cases. Then in the next section, we 
shall complete our analysis by showing that while the geometric picture only makes sense for integer 
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values of m and N, certain formulas and results that one derives using the geometric picture retain 
their validity as m and N are allowed to vary continuously. 

As we have noted, the operators Ka and ICa arose in the study of the Kac model, where they 
measured correlations. We start by explaining the geoemtry behind the simple operator Ka, for 
which the geometric vales of 7 = (A^ — 2)/2, N a positive integer. 



2.1 The geometric origins of the correlation operator Ka 

As is well known, when 7 = (A^ — 2)/2, is simply the image of the uniform probability measure 
(7^ on S^^^, the unit sphere in R^, under the map x 1— t- x • e , where e is any unit vector in M^. 
That is, if eis any unit vector in M^, and / is any bounded measurable function on [—1, 1], then 

/ f{x.e)da^= /(t)d/.«^-')/^)(t) . 
Js^-i J-i 

Let ui and n2 be any two unit vectors in M^, and define a bilinear form on 

by 

qu^,u2if^9)= f{x-ui)g{x-U2)d(TN ■ (2.1) 



We claim that Qui,u2(/)5') is symmetric in / and 5, and depends on the choice of and U2 only 
through a ■.= ui-U2 

To see this, let T be the refiection in about the hyperplane orthogonal to U2 — ui- Then 
T{u2) = ui and T{ui) = U2, and hence, by the invariance of daN under orthogonal transformation 
of R^, 

Qui,u2if,9) = Qu2,uM^9) = QuuuiiaJ) ■ 

A similar argument using a rotation that fixes, say, U2 shows that this bilinear form depends 
on U2 and ui only through a := ui ■ U2, and this established the claim. 

We may now use the quadratic form qu^^^ifid) define an operator Ka where a = ui- U2- It 
turns out that this operator associated is exactly the operator Ka defined in (jl.l6p : 

2.1 PROPOSITION. For any N > 1 and any —1 < a < 1, and all continuous functions f and 
g on [-1,1], 

{Kaf,g)L2,MN-2)/2)^= I f{x-Ui)g{x-U2)d(7N ■ (2.2) 

where Ka is the operator on L'^{fi^"''^) for 7 = (A^ — 2)/2 defined in \1.16\) . 

Proof: This is a calculation based on he following system of coordinates on S^^^: Define 

: 5^-2 X [-1,1] ^5^-1 

by 0(2/,t) = {Vl^yi...,VT — t'^yN-i,t). Evidently for any y G ^ and any t G [—1,1], 
4>{y,'t) S S^~^. It is then easy to check, as in |8], that for any function h on S^~^, 



h{x)d(TM{x) - 
5^-1 J-\ 



H<Piy,i))dcrN-iiy) 

SN-2 



d^((^-3)/2)(^) 
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We now apply this to the integral in (|2.2p . Let {ei, . . . ,e]y} be the standard orthonormal basis for 
M^. Take 

ui = ejv and U2 = clcn + \/l — a'^CN-i ■ (2-3) 
Then with h{x) = f{x ■ ui)g{x ■ U2), we obtain (|2.2p . ■ 

It is clear from (j2.2p that -fCa is self adjoint for 7 = (A^ — 2)/2. Once one knows the self- 
adjointness for these special values of 7, it is natural to seek a direct proof - without lifting the 
functions onto spheres. What one finds is the "intrinsic" quadratic form representation (jl.lSp that 
we gave in the introduction when we proved Gegenbauer's formula. While (jl.lSp may be less elegant 
than (j2.2p . it has the advantage that it is valid for all 7 > 1/2. 

We close this subsection by giving simple probabilistic interpretation Ka which explain our use 
of the term "correlation": Think of S^~^, equipped with dcr^v as a probability space, and think 
of /(x • ui) as a random variable on this probability space. Then, the conditional expectation of 
f{x ■ ui) given 2; • n2 is the function h{x ■ U2) such that E [h{x • U2)g{x ■ U2)\ = E [f{x ■ ui)g{x ■ U2)\ 
for all continuous bounded functions g. Since 

(^a/,5>L2{^{{Jv-2)/2)) ='Ej[Kaf{x-e)g{x-e)] 

for any unit vector e, in particular for e = U2, we see from (|2.2p that 

Kaf{t) = E{/(x • 62) I X • ei = t } . (2.4) 



2.2 The geometric origins of the correlation operator /Cq 

The Jacobi polynomial version of Proposition 12. H leading to /Cq instead of Ka^ is only slightly more 
complicated than the original. Note that any vector x S R'"'^ can be written as an A^-tuple of 
vectors in M"^, x = (xi, . . . , x^r), and hence may be identified with the m x N matrix 

[x] = [xi, ...,xn] (2.5) 

whose jth column is xj. Then for any vector u G M^, the matrix product [x]u is well defined in 
R™. It is easy to see that if x G 5™^"^ and u G S^'^, then [x]u lies in B, the unit bah in M™. 
Therefore, given two unit vectors ui and IL2 in M.^ , and any two functions / and g on B, define 

Qui,u2if,9)= f{[x]ui)g{[x]u2)damN ■ (2.6) 

As before, this will depend only on the choices of ui and U2 through a = ui ■ U2- Hence we may 
use this bilinear form to define a family of self-adjoint Markov operators on L^(dz^m,7v)- Our next 
proposition says that the operators we obtain this way are exactly the }Ca' 

2.2 PROPOSITION. For any N > 2 and m>l, and any -1< a < 1, and all f,g£ C{B), 



{f,^ag)L^{B,u^t,)= f{[x\ui)g{[x\u2)damN , (2.7) 

JgrnN-l 



where ICa is the operator defined in ^\ 
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Proof: We proceed exactly as in the proof of (|2.2p . Define 



by (f){y,v) = — \v\^yi . . . , \/l — \v\'^yN-i,v). It is then easy to check, as in [8], that for any 
function h on S"^^~^, 



h{x)damN{x) 

gl7lN—l 



h{(t>{y,v))darn(N-i){y) 



dl'm,Niv) ■ 



We now apply this to the integral in (|2.6|) with ui and U2 given by (|2.3p . With h{x) = 
f{[x]ui)g{[x]u2), we obtain (fT^ . ■ 
As before, each ICa is a self adjoint Markov operator on L^(dfm,Af)) and has an interpretation 
as a conditional expectation operator: ICa, acting on functions on B, such that for all v £ B, 

fCagiv) = E{g{[x]e2) \ [x]ei = v} . (2.8) 

In the next subsection, we exploit the self-adjointness of ICa to obtain the the product formula 
for Jacobi polynomials in the geometric cases. 

2.3 Spectral analysis of ICa and a product formula in the geometric cases 

In this subsection we study the operator ICa restricted to various invariant subspaces. As we have 
seen, the restriction of ICa to the subspace of rotationally invariant subspaces gives Gasper's kernel 
Kafl ■ The study of ICa on other invariant subspaces leads to the Theorems 11.101 and 11.111 

2.3 LEMMA. For all a G (—1, 1), and all m > 1, N > 2, ICa has the following properties: 

(1) ICa is self adjoint on L'^{i^m,N)- 

(2) If f is a polynomial of degree n on B, then so is fCaf ■ 

(3) For any continuous function f , and any unit vector e, lim/Ca/(te) = fiae). 

(4) For any rotation R on M", /Ca(/ o R) = {Kaf) o R. 

Proof: We argue very much as we did in the ultraspherical case, except of course for the proof of 
(4), which is a new multidimensional feature. 

Propisition 12.21 which expresses ICa in terms of a quadratic form immediately yields (1). As 
for (2), note that di^m,,Af-i(s) is even in s. Therefore, if m is any integer, all of the terms that are 
of odd degree in s that one obtains upon expansion of ^at + sVl — a^\/l — drop out of the 
integral. Hence, what remains is a polynomial in t of degree m. 

Further, (3) follows by the dominated convergence formula; take the limit under the integral 
sign, and use 

lim f (ate + s Vl - aVl - t^) = f{ae) . 

This is independent of s, and since Vm^N-i is a probability measure, (3) now follows. Finally, (4) 
follows from the rotational invariance of Vm,N-i- B 
Since ICa commutes with rotations we can study its action on the irreducible subspaces of the 
rotation group. We begin by considering the action of ICa on the radial functions on B, and shall 
deduce an elegant product formula for Jacobi polynomials directly from Lemma |2.3[ 
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Note that if g is a polynomial of degree at most n in one real variable, and the function / on B 
is defined by f{v) = then by parts (2) and (4) of Lemma 12.31 lCaf{v) is again of this same 

form ~ a polynomial of degree at most n in \v\'^. Thus, the subspace of such functions is invariant 
under fCa- 

Since by part (1) of Lemma I2.3| /Cq is self-adjoint on L^(fm jv)i it may be diagonalized on each 
invariant subspace. It easily follows from here that for each integer n > 0, there is a polynomial qn 
such that with fn{v) = qn{\v\'^), fn is an eigenvector of ICa with eigenvalue An (a), and that the fn, 
appropriately normalized constitute a unit orthonormal basis for the subspace of radial functions 
in L'^{i'm j\i). By the explicit form of i^m,N given in (|1.24p . 



where Cm,N is a normalization constant. Making the change of variables t = r"^, one now recognizes 
the Qn as being the Jacobi polynomials for a, /3 given by (|1.26|) . translated and scaled so the domain 
is [0, 1] instead of [—1, 1]. 

To determine the eigenvalues A„(a), apply the evaluation formula, part (3) oi Lemma 12.31 to 
see that \nia)qn{l) = lim|^|^i A„(a)gn(bP) = fCafniv) = Qnia'^)- That is: 



Since the restriction of ICa to the radial functions on B is clearly positivity preserving, and clearly 
preserves the constants, it follows that the A„,(a) := '?n(a^)/Qn(l) are a Markov sequence for the 



Thus, condition (1) of Theorem 11.51 is satisfied, and as a consequence of Theorem 11.51 we 
have therefore solved the Markov sequence problem for the Jacobi polynomials, and have proved 
a product formula for them, in the geometric cases. One can of course undo the scaling and 
translation, and write this all out explicitly for the usual Jacobi polynomial defined on [—1, 1]. The 
resultis, of course. Gasper's product formula. We shall do this, but first notice that there is more 
to be obtained from the analysis of ICa- So far, we have only considered the restriction of ICa to 
the radial functions. The spectral analysis of ICa on other invariant subspaces prvides additional 
formulas identifying ratios of Jacobi polynomials as eigenvalues of self adjoint operators. We shall 
use these formula (and their extension to general values of a and P) to prove Theorems 11.101 and 



2.4 The spectral analysis of ICa on non-radial functions 

For each integer £ > 0, let Hi denote the space of harmonic polynomials on that are homoge- 
neous of degree £. Restricted to B, the functions in constitute a closed subspace in L'^{i'm,N), 
which we again denote by 

For each £, Tii is an eigenspace of ICa- In fact, for each H G H^, 



That is, the restriction of ICa to i7 E Hi is times the identity. One way to see this is to use the 
mean value property of harmonic functions and the formula (jl.24p . Since the measure di'm,N-i is 
radially symmetric, we see that ICaH{v) = H{av), which, by the homogeneity, is a^H{v). 





[mi 



ICaH{v) = a^H{v) . 



(2.9) 
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There is another more algebraic argument that tells us somewhat more: 

2.4 LEMMA. The spectrum of K-a is discrete, and its eigenf unctions are of the form g{\v\'^)H(v), 
where g{\v\'^) is a polynomial in \v\'^ and H G Tii for some £. Moreover, if g{\v\'^)H(v) is an 
eigenfunction, then, so is g{\v\'^)H{v), for any non zero H £ Hi. 

Proof: By Lemma 12.31 the operator fCa leaves the space of polynomials of degree n invariant for 
any n. Hence, by the Weierstrass theorem the eigenfunctions consist of polynomials. Further, since 
fCa commutes with rotations, any eigenfunction must be of the form 

F{v) = f{\v\)y'{£^ (2.10) 

where is a spherical harmonic, i.e., (|~~[) ~ \^\~^Hi[v) where Hi{v) is a homogeneous har- 
monic polynomial of degree I. We have to show that is a polynomial in v, i.e., a 
polynomial of the variable \v\'^ . 

n 

Since F{v) is a polynomial of degree n we can write it as F{v) = q-miv) where qm{v) is 

m=0 

homogeneous of degree m. In turn, each of these polynomials can be expanded in terms of homoge- 
neous harmonic polynomials, i.e., qm{v) = H^iv) + \v\^Hm-2{v) + \v\'^Hm-i{v) + • • • ■ This shows 
that 

n 

F{v) = Y,9k{\v\^)Hk{v) (2.11) 
fc=o 

for some polynomials gk- The result follows from (I2.10p and the orthogonality properties of the 
spherical harmonics. The final statement follows from Schur's Lemma since tCa commutes with 
rotations and rotations act irreducibly on T-L^. ■ 
Now, since polynomials on [0,1] are uniformly dense in C([0, 1]), it follows from the Lemma 
(and the fact that fCa is Markov) that for each function g G C([0, 1]), and each H G T-Li, and all 
a G (—1, 1), there is a G C([0, 1]) so that 

}Caf{v)=ga{\yf)H{v) where f{v)=g{\v\^)H{v) (2.12) 

with H G Tie being the same on both sides. The transformation g >-^g is clearly linear, and as one 
sees from the proof of Lemma [2.4| independent of the choice of H. We now use this transformation 
to generalize the definition of the operator in ()1.2ip . 

To make efficient contact with the theory of Jacobi polynomials, it is better to write our radial 
functions in the form /i(2|up — 1) instead of ^(It^p). For any non zero H in any H^, we define Vh 
to be the subspace of L?'{um,N-i) consisting of functions of the form 

f{v) = h{2\v\^-l)H{v) , 

where is a function on [—1,1]. We then generalize the definition (|1.21|) as follows: 
For each £ > 0, fix some non zero H G T-Li- Then for h G C([— 1, 1]), define Ka/h by 

{Ka,ih){2\v\^ - l)Hiv) = {JCafKv) where fiv) := hi2\v\^ - l)H{v) . (2.13) 
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By the last statement in Lemma 12.41 [Ka/h){2\v\'^ — l)H{v) does not depend on the particular 
choice of H in Hg. Further, by Lemma 12.41 the eigenfunctions of fCa are of the form 

where hn^t is a polynomial of degree n. The eigenfunctions of Ka/ are the polynomials We 
now identify these as Jacobi polynomials: 

First, we normalize our choice of H G He so that / \H{v)\'^dam = 1- Then for any positive 

integers n ^ p, integrating in polar coordinates and using (|1.24p we find 

= {fn/,fp,i)LHu^j,) = [ VK2kl' - l)V(2|t>|2 - l)\H{v)\^dl^^,Niv) 

JB 

|^m(Af-l)-l| rl 



-1| 



V(25' - 1) V(2s' - 1)(1 - 52)MiV-l)-2)/2^2^+™-l^^ _ 



Making the now familiar change of variables t = 2s^ — 1, we find 

-1 

which is the orthogonality relation defining the Jacobi polynomials pl^'^~^^^ with a and /3 given by 

We now determine the eigenvalues Xn,eia) such that Ka/hn/ = \n/{o)hn,t- By (I2.13p . if we 
define / by f{v) = hn,e{2\v\'^ — l)H{v), we have Xn/{a)f{v) = JCaf{v). Then, for any unit vector e 
in M*" with H{e) ^ 0, we have from part (4) of Lemma [2.3l that Xn/{a)f{e) = f{ae). which, by the 
homogeneity of H, means that Xn,e{a) = a^hn/{2a'^ — l)/hn/{l)- We summarize our conclusions in 
a lemma: 

2.5 LEMMA. Fix dimensions m > 1 and N > 2, and let a and f3 he given by \1.26\) . Then for 
each integer i > and each a € [—1, 1], the operator Ka/ is self adjoint on L'^{fi^'^'^'^^^), and is 
diagonalized by the Jacobi polynomial basis {p^'^~^^^}n>o- Moreover, the corresponding sequence of 
eigenvalues {A„/(a)}n>o is given by 

^nA-) = ■ (2.14) 

Pn \'-) 

2.6 REMARK. It is worth remarking that the operators Ka/ are not positivity preserving for 
£ > 0. Nonetheless, the eigenvalues Xn,i{a) are eigenvalues of a Markov operator, namely ICa- 

3 General values of a and ^ 

Our next goal is to extend this analysis of the previous section to arbitrary values of a > /3 > —1/2. 
We seek a direct expression of Ka/, not explicitly involving ICa, so that we may then freely vary 
the dimensions. The following notation will be useful: For v and y in B, define 



w{v, y, a) = av + — a^\J\ — \v^y 
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Then by Lemma 12.21 for any / of the form f{v) = h(2\v\'^ — l)H[v), 

K^af{v)= h{2\w{v,y,a)\'^ -l)H{w{v,y,a))di^rn,N-i{y) ■ (3.1) 

JB 

To proceed, we now make a judicious choice of H to simplify the computations before us: 
• We choose H so that H{v/\v\) is the zonal spherical harmonic of degree i with the axis along the 
unit vector e in M™. That is, 

H{v) = \v\^p\^"''^'^''^\e- v/\v\) = \v\^pf\e- v/\v\) , (3.2) 

where, as before, p^j^^ is an ultraspherical polynomial, and /3 is again given by (jl.26p . The advantage 
of our particular choice of H is that H{w{v, y, a)) depends on w{v, y, a) only through {w{v, y, a)) ■ e 
and \w{v,y,a)\'^ . Specifically, 

H{w{v,y,a)) = \w{v,y,a)fpf'' f iwiv^V^a)) e \ ^ 

V \w{v,y,a)\ J 

Note that 

w{v,y,a) -e = as + \/l — a^\/l — s^r cos 9 (3-4) 

and 

\w{v,y,a)\'^ = a^s^ + (1 - a^)(l - s^y"^ + 2a\/l - aVl - s'^rscosO (3.5) 

where s = r = |y|, and v ■ y = srcosO, so that ()3.2p can be written as an integral over r and 6, 
using the measure defined in (|1.22p . These are the coordinates we used in the proof of Lemma 11.91 
to "liberate" the values of a and /3 in the i = case. 
By (j2.13p and (j3.2p evaluated at v = se, we have 

{Ka,ih) (2s2 - l)s^H{e) = [ hi2\w{se,y,a)\^ - l)H{w{se,y,a))di^^,N-iiy) ■ 

JB 

For the particular choice of H made in (|3.2p . this reduces to 

{Ka,ih) {2s' - 1) = S-' [ h{2\w{se,y,a)f - l)\w{se,y,a)fpi^^ / (u;(^e,j/, a)) • e \ ^ 

Jb \ \w{se,y,a)\ ) 

(3.6) 

since p^P {x) / p^P {!) = P^^\x), the ultraspherical polynomial normalized by the condition 
Pr(l) = l. 

Next, since the integrand depends on only on s, r and cos 9, we can use (ll.24p and (I1.22P to 
write this in terms of an integration against drria^p with a and (3 related to m and through 

{K,^,h) {2s'-l) = S-' r rh{2\w{se,y,a)\'-l)\w{se,y,a)fpi^^ A^ge,j/, a)) • e \ ^ ^ 
Jo Jo V \w{se,y,a)\ J 

(3.7) 

3.1 DEFINITION. For aU a > /3 > -1/2, we define Ka/ by formula ([ST]). By the calculation 
just made, this coincides with the definition made in (j2.13p for a and (3 satisfying (jl.26p . 

The next Lemma gives a more explicit formula for Ka/- 
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3.2 LEMMA. For all a > (3 > —1/2, and all non negative integers i, 



Ka/h{t) 




h 



^0 

M 

j=0 \ 



a2(l + t) + 62(1 _ t)r2 + 2abVl - t^r cos d 



a'^^ibry 



/ 1 - t 



i/2 



1 +t / J ^ 



dmQ,,/3(r,6') , 



(3.8) 



where b = \/l — as before. 



Proof: First make the change of variables t = 2s^ — 1 in the defining formula (13. 7p . Under this 
change of variable, (13.41) and (13.51) become 



and 



, y,aj ■ e = ayl — h oa/ — - — f cos d 



.l+t ,nl-t 



\w{v, y, a)|2 = — 1- 6^ — - — + abyl — l^r cos 9 , 



(3.9) 
(3.10) 



and we deduce from (|3.7|) that 



l + t 

2 



h 



JO 



^(1 + t) + 6^(1 - t)r'^ + 2ab\/ 1- f^r cos e 



a-(l + i) + 6^(1 - t)r2 + 2a6\/l - cos 6* 



^/2 



X P. 



+ 6V1 - tr 



cos I 



x2 (1 + 1) + 62(1 _ i)^2 _^ 2a6Vl - cos 6* 



Y/^ dm„,;3(r,6') 



(3.11) 



The Laplace formula for the ultraspherical polynomials [21j, p. 94, which is a simple consequence 
of Gegenbauer's identity, can be written as 



P<«(„_r(/i + i/2) r 



V5T(/3) 



j [x + — 1 



COS sm 



-,2/3-1 



(3.12) 
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With X = - — J- — ^ r-j- , we obtain, using the binomial formula, 

\w{v,y,a)\ 



1 +t) +6^(1 - t)r^ + 2abV I- f^r cose 



X I, aVTTt + bVT^tr cose _| 



\a 



'l + t) + 62(1 - t)r2 + 2a6Vl - cos ( 



{aVTTtY-^ibrVT^y (3.14) 

i=o V / 

^^^wYx^^ /"(cos 6* + \/cos2 - 1 cos sin^^-i (3.15) 
V^r(/i^) Jq 

^ ( M (a7rTi)^~-'(6r^/r^)^pf )(cos0) (3.16) 



j=0 



Proof of Theorem ll.lOt By choosing h = Gasper's product formula (Theorem ll.lOp 

follows immediately from Lemma l3.2i ■ 
Our next goal is to show that for all a > /3 > —1/2, Ka/ possesses the crucial properties of 
self-adjointness, polynomial preservation and the limiting value identity that it inherits from /C^ 
when a and /3 satisfy (|1.26p . 

3.3 LEMMA. For all a G (—1,1), a > /3 > —1/2, and integers £ > 0, the operator Ka/ on 
C([— 1,1]) as defined in \1.21\) has the following properties: 

(1) Ka/ is self adjoint on L'^{ii°^'^). 

(2) The space of polynomials of any fixed degree is invariant under Ka/. 

(3) For any continuous function h, \im Ka eh(t) = a^h{2a^ — 1). 

i— 5-1 ' 

Proof: It is obvious from (13. lip that limt_j>i Ka/h{t) = a^h{2a^ — 1), property (3) is taken care of. 
Next consider the polynomial preservation, property (^^j. It suffices to show that for each natural 
number n, if h{t) = [t + 1)", then Ka/h{t) is a polynomial of order n. 
For this choice of h, 



h p(l + t) + 6^(1 - t)r'^ + 2ah\J\ - t^r cos 6* 



1 



y -, {a^{\ + t) + 6^(1 - t)Y"" ( 2a6\/l-tVcos(e) 

^ (n- mV.ml v ^ ' ^ V 

m=l ^ ' 

Thus, from Lemma 13.2^ Ka/h{t) is a sum of multiples of terms of the form 

(-1 _ f\ pi p-IT 

1 + Jo Jo ' 

where Q{t) is a polynomial of degree n — m. Then, be the orthogonality properties of the ultras- 
pherical polynomials, 

/•TT 

cos"" epl^'^ (cos e) sin^^ ede = 
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unless m + A; is even and m> k, in which case (1 — t"^)^/"^ I | = — — — - is a polyno- 

- ' ^ ^ Vl + V (1 + 1)("^-'=)/2 ^ ^ 

mial of degree m. Thus, for this choice of h, Ka/h{t) is a sum of terms each of which is a polynomial 
of degree n, and thus (2) is proved. 

We next deal with self-adjointness. To see this in a simple way, we do not use the formula for 
Ka/ given in Lemma 13.21 but instead work directly from the expression p.7p . We shall show that 
the bilinear form 

q{hi,h2) := 2c,,;3 [' hi{2s^ - 1) {Ka/h2) {2s^ - 1)(1 - (3.17) 
JO 

is symmetric. This is easily seen in case a and /3 related to m and N through (|1.26p since then 
with fj{v) = hj{2\v\^ — l)H{v), j = 1,2, easy computations reveal that the right hand side is a 
constant multiple of {fi-,K,af2) l'^{v^ ^) To see this in general, we proceed exactly as in the proof of 
Lemma 11.91 making the same sequences of coordinate changes 

{r,e) {x,y) {x ,y) {p,(j)) . 

Under this sequence of changes of variables, \w{v,y,a)\ becomes simply p, as we have seen in the 
proof of Lemma [L9l and w{v,y,a) ■ e/\w{v,y,a)\ becomes simply cos{(p), as simple computations 
reveal. Then, with q{hi,h2) defined in (j3.17p . we find that 

q{hiM) = 2c„,/36-2a /■\^(2s2 - l)/i2(2p2 _ 1) 

JO JO JO 

X (62 _ ^2 - p2 + 2aps COS 0)"-^" ^ (cos sin^/^ 0d</) dpds , 



(3.18) 



This takes care of property {1). 

With this lemma in hand, we now easily extend Lemma 



3.4 THEOREM. For all a > P > -1/2, all integers i > 0, and each a G (-1,1), the operator 
Kfi/ is self adjoint on L2(^{"./3+^))^ qj^^^ j^g diagonalized by the Jacobi polynomial basis {pn''^^^^}n>o- 
Moreover, the corresponding sequence of eigenvalues {A„^^(a)}„>o is given by 

Pn 



Proof: This is an immediate consequence of Lemma 13.31 and Theorem 11.81 
Proof of Theorem II. lit By Theorem 13.41 



pt'^'\ t) 



%t^+'H2a' - 1) = a-' (K^,ept^+'A (t) . (3.20) 
1) ^ ^ 



Now, the left hand side is a polynomial in a, and so is the right hand side. Hence we may 
extend the range of a from [—1,1] to all of M. Since in (j3.8p . b stands for \/l — a2. All odd terms in 
Vl — must drop out of when the integration is made, and for a > 1, we will get the signs right 
if we replace b = \/\ — with i\/ — 1. 
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Doing this, and then dividing both sides of (j3.20p by (2a'^)", and taking the hmit a — )• oo, only 
the leading terms in the Jacobi polynomials contribute, and we obtain we obtain 



Fn 



it) 



1 

^0 
X 



(l + t)- (1 -t)r 



E 

.fc=0 



+ i\/l — t^r cos 



l + t 



ir)'=pf^(cosl 



dma,[s{r,9) 



This is the formula in Theorem 11.111 



4 Bounds on the extremal eigenvalues and convergence of 
associated eigenfuncton expansions 

Our objective in this section is to obtain bounds on the magnitudes of the eigenvalues in the 
extremal Markov sequences that govern the way these magnitudes decrease to zero as n increases. 
We start with the case of the ultraspherical polynomials. 

4.1 THEOREM. For all 7 > 0, and all -1 < a < 1, 

2c, 



< 



1 - a2)7 1 2/ 



(4.1) 



Moreover, if p > I/7, {KaY, the pth power of Ka, is trace class. This criterion for belonging to the 
trace class is sharp in that for a = 0, where exact calculuations are simple, one finds Kq is trace 
class if and only if p > I/7. 

We shall prove an analog of this Theorem for Jacobi polynomials, and the proof will be quite 
similar. Therefore, before plunging into the details, we explain the strategy. 
The starting point is the Laplace identity (13.12P which can be written as 



pl^\a) 
i^^7^(l) 



— = / (a + isy/l 
1) J-i ^ 



(7-1/2) 



is) 



Observe that 



< 1 . 



|o + isVl - a2|2 = 1 - (1 - a2)(l - 

Fixing a, define Cx to be the subset of [—1, 1] on which \a + isVT" 
from (j3.12p and layer-cake that 



i2|2 > 1 



(4.2) 
A. It follows 



P^7^(a) 



pi''^(l) 



< 



J a + is^/T 



(7-1/2), 



s) < 



n 



(1 - A)("-2)/2^(7-i/2)(c';,)j^ _ ^4 3) 



Hence, an estimate on the rate that fi^'^~^^'^\Cx) decreases to zero as A decreases to zero yields 
a bound on the rate at which \pn\a) / p!^\l)\ decreases as n increases. This will yield us bounds 
that hold uniformly in a in any compact subset of (—1,1). While we are ignoring phase cancelations 
in the estimate (j4.3p . there are no phase cancelations for a = 0, and an exact calculation gives the 
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same n decay. Thus no better bound can hold uniformly in a on closed symmetric intervals of 
(-1,1). 

We prove a bound on ii^'^~^/'^\Cx) in the next lemma, and then proceed with the proof of the 
theorem. 



4.2 LEMMA. n'-'^'^^'^\Cx) < 2c^_i 



Proof: Note that from g2D, + > 1- A {l-t'^){l- s^) < \. Hence, for s € Ca, 

1 - < A/(l - t^), and therefore, 



(7-1/2) 



(Ca) = 2c^„i/2 [' (1 - s'r-'ds < 2c^_i/2 (rr^) 



Proof of Theorem I4.lt Applying Lemma 1^2] in (j4.3p . we obtain 



P^7^(i) 







c^_i/2 r(§)r(7 + i) 

-n- 



(i-t2)7 r(§ + 7 + i) 



(4.4) 



Then since (1 — e *)''' < s'^ for 7 > 0, 

r(5) 1 



r(f+7 + i) r(7 + i) 



°° n /'n\-(7+l) 



(4.5) 



Combining this with (j4.4p we obtain the bound ()4.ip . 

Finally, we consider the case a = 0. Then there is no phase cancelation, and one readily 
computes 



p£^(i) 



c^_i/2r(n + ^)r(7) 

V2 r(n + 7+i) 



~ n 



and we see that the n~'^ bound on the absolute value of the extremal eigenvalues is the best possible 
that can hold uniformly for o in closed intervals of (—1, 1). ■ 
We now turn to the analog of Theoem 14.11 for Jacobi polynomials. 



4.3 THEOREM. For all a > (3 > -1/2, all a £ (-1, 1), all n > 0, and all £>0, 



Pn 



(a) 



Pn 



(1) 



< 



1 + 



1 - a 
1 + a 



1/2 



K„,^(a)r(a + -j (- 



(4.6) 



where 



Thus, for p > 1/(q + 1/2) > 0, {Ka/Y is trace class. 

4.4 REMARK. The exponent on n is determined by a alone; it is independent of fi and i. 
We begin with a lemma that is the analog of Lemma 14. 2t 



(4.7) 
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4.5 LEMMA. For each fixed a G (-1,1) and A > 0, define C\ to he the subset of [0,1] x [0, vr] 
given by 

Cx ■.= {{r,9) : R'irJ) > 1-A } . (4.8) 



where 



R{r,t 



(1 +a) - (1 - a)r2 



a^r cos ( 



(4.9) 



Then the measure of Cx with respect to dnia^p satisfies 

m^,p{Cx) < i^«,;3(a)A°+i/2 . 
where Ka^pia) is given by j^. 7] ). 

Proof: Define A = (1 + a)/2 and B = [1 — a)/2. Then we can write 

R^ir, e) = A^ + B'^r^ + lABr"^ cos(26') . 

Since A + B = 1, 1 - R^{r, 6) = 5^(1 - r^) + 2AB{1 - r^ cos(26')). This is a sum of positive terms, 
and so for any A > 0, whenever 1 — i?^(r, 6) < A, we have both 



B^{l-r^)<X and 2^5(1 - cos(20)) < A 



(4.10) 



The first of these conditions imphesl— < A/(-B^(l+r^)(< A/i3^, and then 1—r < (1 — r)(l+r) 
1 — < X/B'^, so that r > 1 — X/B"^. Thus, every where on Cx, 



1-r'^ < — 
^2 



and 



r > 1 

B2 



Next, we turn to the second condition in (j4.10p . This can be written as cos(20) > 1 — X/{2AB), 
which certainly imphes cos (20) > 1 — X/{2AB), which imphes that sin^ 6 < X/{AAB). FinaUy, since 
on [0,7r/2], {2/tt)0 < sm{9), with a similar estimate on [7r/2,7r], the second condition in ()4.10p 
implies that either < 6* < {tt / 4) y^X/{AB) , or else vr - {tt / 4) ^/X/{AB) < 6^ < vr. 



Altogether then, {(r, 61) : i?2(r, 6*) > 1 - A } is contained in [1 - A/S^ 1] x [0, (7r/4)yA/(^5)] U 
[1 — X/B'^, 1] X [vr — (7r/4)y A/(Ai?), vr], and moreover, everywhere on this set, 

A 



2 A 
1 -r^ < 



and 



sin^ d < 



" - 

Integrating over the two rectangles using the above bounds yields the estimate 

mo^ACx) < C„,/3Vr-l23-2/3i?-(2-/3+l/2)^-(/3+l/2);^a+l/2 _ 

Replacing A and B by their definition in terms of t, one obtains the bound (|4.9p . 

Proof of Theorem 14. 3t The starting point is Theorem 11.111 which provides the identity 



Jo 



[l + a] 



[1 — a)r 



+ i 



a^r cos ( 



E 

.fc=0 



. xk/2 

i^j (^r)'^pf (COS0) 



dmo,,i3ir,B) . 
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By the definition of R{r, 9) in Lemma |4.5| and the fact that Pn\x) < 1, we obtain 



< 



1 rn 



Jo 

1 PTT 



< 



Jo 



1 + 



1 + 



Rir,9y 

R{r,ey 

1 - a 
1 + a 

1 \ 1/2' 
1 — a ^ ' 



E 

.fc=0 



fc/2 



/c / \l + aj 



dma,i3{r,e) 



1/2- 



1 /-TT 



1 + a 



JO 
1 







Now applying Lemma I4.5[ and then estimating the ratio of Gamma functions as in 



Pn 



(a) 



< 



< 



1 + 



1 + 



1 + 



1-a 



1/2 



1 + a, 

1 \ 1/2' 

1 — a ^ ' 



1 + a 
1 — a ^ ' 



^K^Aa) / (1 - A)('^-2)/2A"+V2dA 

JO 

-Aa,/3(aj 



r(f + a + |: 



1 + a 



. / 3\ /n\-("+i/2) 
i^a,/3(a)r(a + -J (- 



(4.11) 



(4.12) 



We may apply these results to study the convergence of the eigenfunction expansions for the 
operators Ka and Ka^e- Indeed, since eigenvalues of Ka^e are 

,pt'^'\2a'-l) 



and the eigenfunctions are the p^'^^^^ , the formal eigenfunction expansion of the Kernel for Ka , 



IS 



e^ pt'^'\2a--l)pt'^'\x)pt^^'\y) 



The eigenvalue bounds obtained above can be used to show that for a G (—1,1) and a > 1/2, 
these formal series actually converge uniformly for x and y in compact intervals of (—1, 1). To do 
this, we need bounds on the eigenfunctions as well as the eigenvalues. Because of the close relation 
between the eigenfunctions and the eigenvalues in this context, we could obtain the eigenfunction 
bouts from Theorem 14.31 but it will be instructive to obtain these instead from a well known but 
deep result of Nevai, Erdelyi, and Magnus. [2DJ: For all a > —1/2 and f3 > —1/2 and all non 
negative integers n. 



max. 



:,el-i,i]Vl-xMx)p'^'Hxf < 



2e(2 + ^/a^Tp^) 



IT 



(4.13) 
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Thus, for each a G (—1, 1) and each r < 1, there is a constant C such that 

uniformly for x,y £ [— r, r]. With the a = /3 cases coming from Theorem 14. 1 1 and ()1.4p . this proves: 
4.6 THEOREM. For all q > 1/2 and a > /3 > -1/2, and all -1 < x,y,z < 1, the sum 

converges absolutely and uniformly on compacts, and the operator whose kernel the sum defines is 
trace class. 

For £ = 0, this is the eigenfunction expansion of Gasper's operator Ka^, which is Markov. Thus, 
for a > 1/2 and a > f3 > —1/2, where the sum in (I4.14p . converges pointwise, it defines a kernel 
that is pointwise positive.. 



5 Historical Remarks 

In this section we give a brief discussion of work done on the Markov sequence problem for Jacobi 
polynomials by Gasper, Koornwinder and Askey with the aim of clarifying the context of the present 
paper. 

When Gasper took up his work on the Markov sequence problem for Jacobi polynomials, the 
main obstacle was the lack of an analog to Gegenbauer's indentity. Therefore, Gasper worked 
backwards towards one: It is clear that the operator Kafi on L^(^("''')) defined by 

K^,oHx) = E (at. . ^ pt'Hx) / pt'\ymyW^-''Hy) (5.i) 

n=0 Pn ' (1) -^-1 

for a G (—1,1) is self adjoint, satisifes Kafll = 1, and has {pi^'^\z) / pl^'^\l)}n>o as its sequence 
of eigenvalues, so that if there is to be an analog of Gegenbauer's identity for Jacobi polynomials, 
it must refer to this operator. 

What is not at all clear from the eigenfunction expansion is whether or not Kafi preserves 
positivity, or whether Kafi even has a kernel Kafi{x,y), which formally would be 



oo (a 



''^\2a^ - l)pt^\x)pt^\y) 



Kafl{x, y) = 2_^ .o) ■ (5.2) 



n=0 



pt'\l) 



Gasper's Theorem as stated above was proved by him in |llj . where he evaluated the sum. This 
is rather involved, but here is a brief sketch: In [TT] Gasper formally defines a kernel G{x, y, z; a, /3) 
as 

oo a, fir r)i\ 

G(cos2</>,cos2^,cos2e;a,/3) = Y,K'^ ^'" f °' /^ ^'^(cos 2V^)K'^(cos 2g), 

n=0 (1) 
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where hn'^ is the square of the inverse of the norm of pn'^ ■ Then using a formula of Watson, [22] 
p. 413, he shows that this sum of triple products of Jacobi polynomials is equal to an integral of a 
triple product of Bessel functions with the restriction a > —1/2,13 > —1/2 and cos6 ^ \ cos (V'i </>)!• 
When a > f3 > —1/2, he was able to evaluate the integral of the triple product of Bessel functions 
with the result that 

G(cos 2(j), cos 2-0, cos 29; a, 13) 

_ r(a + l)(sin0sin?/;sin(9)-2« ^ 2 



„x^,^ ,x I (1 — cos 6 — COS ■0 — cos 6+ 

2a+/3+ir(a_/3)r(/3 + l/2)r(l/2) 7o 

2 COS cos ip cos 6* cos 7)""'^"-^ sin^^ 7^7 (5.3) 

(COS^ (j) ~\- COS'^ Ip ~j~ cos*^ — 1 \ 
I , or TT depending on whether sin^ (p sin^ ip is less 
2 cos (j) cos ip cos J 

than between, or greater than the two numbers (cos cos '0 i cos^)^. from this Gasper concludes 



that G is non- negative. Then in [12], using the evaluation of the triple integral in terms of hy- 
pergeometric functions. Gasper is able to show that G is nonnegative if a > /3 > —l,a > —1/2, 
and either /3 > —1/2 or a + /3 > 0. Comparing this with equation (jl.lSp we see that the kernel 
is positive so that Gasper's result lays the foundation for a convolution structure associated with 
Jacobi polynomials. 

Later Koornwinder [l^ gave another proof of Gasper's Theorem in the case a > f3 > —1/2. 
Here Koornwinder defines the kernel G as the integral (|1.2ip given above, then he uses his Laplace 
type integal representation for Jacobi polynomials and duality to show that the kernel is equal to 
the triple sum of Jacobi polynomials. Koornwinder obtained his Laplace type formula using group 
theoretic methods and Askey [1] gave a simple analytic proof using Bateman's integral relation 
between hypergeometric functions. The fact that the kernel is continuous and of bounded variation 
allows Koornwinder to show, using the Dirichlet-Jordan test |23j p. 57 and the equiconvergence 
of Jacobi series and cosine series |21] p. 246, that for a > (3 > —1/2 the triple sum converges 
uniformly on compact subsets of < 0, 0, < ^. Later Koornwinder and Schwartz [18] extended 
these results polynomials orthogonal on the biangle, triangle, and simplex. 

In contrast our approach is much more in the spirit of Bakry and Mazet |3], in which they 
solved the Markov sequence problem for various systems of orthogonal polynomials by applying 
functional analytic techniques to certain well chosen kernels of self-adjoint operators. Likewise, our 
approach starts with the construction of the family of operators the operators tCa- The motivation 
for considering the family /C^ comes from previous work on the Kac model [H], [S]. In particular, for 
the restricted parameter values discussed in Section[2l K-a is easily seen to be self adjoint, to preserve 
polynomials and to enjoy the evaluation property, and hence, by Theorem 11.81 its eigenvalues can 
be expressed as ratios of Jacobi polynomials. Restricting this operator to the radial functions gives 
us the operator Kafi, at least for the half integral values of a and (3 given in (jl.26p . In a further 
step we extend the operators i^a,o to the full range, and obtain Gasper's kernel. i^a,o is shown to 
be an extremal Markov operator from which Gasper's product formula follows 

Moreover, the same can be done with the operators Ka^t that appear as restrictions of the 
operators tCa to the other invariant angular momentum subspaces, and in this way we obtain 
Theorems fTTOl and fTTTl 

One final remark is that the convolution structure for Jacobi polynomials can be used to show 
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that Jacobi polynomials form a strong polynomial hypergroup (see Bloom and Heyer [5] for a 
definition of hypergroups). A sequence of orthonormal polynomials {/«} with /o = 1 such that 
11.71 holds is said to have the hypergroup property (see Bakry and Huet [3j). A result of Connett 
and Schwartz [10] (see also Theorem 3.4 in |5j) essentially shows that the only unit orthogonal 
polynomials systems in one dimension that satisfy 11.71 are, up to a translation and scaling, the 
Jacobi polynomials. This can be proved by showing that if a sequence of orthonormal polynomials 
satisfy 11.71 then after translation and scaling they also satisfy the differential equation satisfied 
by the Jacobi polynomials. The argument does not, however, provide another proof of Gasper's 
Thoerem asserting that the Jacobi polynomials do indeed satisfy a product formula. 

A The parabolic biangle and triangle polynomial product formu- 
las 

The parabolic biangle is the region 

B = {(xi,X2) -.0 <xl<xi <1} , 
and the polynomials rn',m{xi,X2) that are orthogonal on B with respect to the measure 

(1 - xi)°(xi - xlf-^/^dxidx2 
can be written, with our conventions as 

C,^ix^,X2) = xfpt^+-H2x, - l)p(f)(^) . (A.l) 

The total degree is n + m. 

The product formula Koornwinder and Schwartz [18] is then given by 



rn\m{xl,X2)rn\m{yl-,y2) 
<d(l,l) 



/ rZ:i{E\EG)dv''^^{T, h,t2,h) (A.2) 



where 

du"''^{r,ti,t2,h) = dfi^{t2)dfi'^{ts)dm^^,3{r,ti) , 

where fi^ is given by ll.lt m^^^ by 11.221 with ti = cos0. Here / = [0, 1], J = [—1, 1]. The symbols E 
and G are given by 



E = {xjyl + (1 - xl){l - yiy + 2xiyi(l - xlf'\l - ylf'^h) 
G = D{C, D{x2/xi,y2/yi; 1, t2); 1, ts) 



1/2 



where 

D{xi,yi;r,ti] 



G 



E{xi,yi;r,ti) 
and, generally, 

D{x, y; r, t) = xy + {I - x^)'^^'^{l - y'^f'^rt . 
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In order to prove this product formula ()A.2p using the method developed here, we must analyze 
the operator 

{JCy,,y,h){xi,X2):= [ h{E\EG)diy'''^{r,tut2,t3) (A.3) 
Ji,J^ 

on the Hilbert space T-L given by the inner product (•, •) defined above. Direct calculations that 
lead to the proof of self-adjointness of the operator K,y^^y2 defined above seem to be very involved. 
A substantial simplification is achieved by writing JCy^^y^ in terms of Gegenbauer operators. For 
the triangle case the Gasper operator will be used instead. It is an immediate consequence of the 
following theorem. 

A.l THEOREM. The operator fCy^^y^ is a self adjoint linear operator on %. It is positivity 
preserving, preserves the function 1 and the space of polynomials in the two variables xi,X2 of a 
given degree. Further we have the evaluation formula 

. ,}^^, , J^yuy2h)ixi,X2) = h{l,l) . (A.4) 

Proof: Clearly, K,y^^y2 is positivity preserving and preserves the function 1. The evaluation formula 
follows by noting that E and G tend to 1 as (xi,X2) tend to (1,1) in the biangle. To see the 
statement concerning polynomial preservation it suffices to prove it for a general monomial x^X2- 
That is we have to show that 

I E^^iEGTdv'^'^irMMM) 
Ji,.P 

is a polynomial in the variables (xi, X2) of total degree less than or equal 2m + n. We shall use the 
fact that the measure d^"'^ is even in ti,t2,t^. Now, 

G" = [CD{x2/xu 2/2/2/1; 1, t2) + (1 - G^fl\l - Z5(x2/xi, 2/2/yi; 1, t2fY'h; 



Y.\l] C''-''D{x2/x^,y2/yrAMr-\l-G^fl\l-D{x2/x^,y2/yr,lM?f'^t. 



k=0 



and integrating this expression with respect to the ^3 we see that only terms with even k contribute 
and we obtain an expression of the form 

[n/2] 

Ycin,2k)C^-'^''Dix2/xuy2/yi;l,t2r-'\l-GY{l-Dix2/xuy2/yi;l,t2f)^ 

k=0 

where c{n, 2k) are positive coefficients. This expression can be rewritten as 

[n/2] k / \ 

J2cin,2k)G^-^\l-GYY.[ {-irD{x2/xi,y2/yi;l,t2r-''^'P . (A.5) 

k=0 p=0 \ ^ / 

Applying the binomial formula to the expression D{x2/xi,y2/yi', 1, ^2)""^'^^^^ and integrating with 
respect to t2 leaves us with a polynomial in the variables {x2/xi)'^ if n is even or is of the form 
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X2/X1 times a polynomial in (x2/xi)^ otherwise. Moreover, it has degree not larger than n. The 
remaining terms are, when multiplied by ^E^m+n^ q£ ^^le form 

which when integrated over ti yields a polynomial in xf if n is even or is of the form xi times a 
polynomial on otherwise. It has degree not larger than 2m + n. Thus, after performing the 
integration over the variables r, ii, t2, is one obtains a sum of terms of the form 

— IJ = x^'^ '^^X2^ n even , with 2q < 2m + n 

or 

— I I = x^^"^^ x^'^^ n even , with 2q + l < 2m + n . 
x^ J 

Thus we obtain a polynomial of the form •p{x\^ X2) whose total degree is not larger than the one we 
started with. It remains to show selfadjointness. It is convenient to write the inner product 

(/,5):= /'dxi(l-xi)° dx2ix^ - xlf-'/''f{xuX2)gixi,X2) . (A.7) 

Jo J-^/xl 

in terms of the functions 



F{p,s) := f{p^,ps) , G{p,s) ■.= g{p^,ps) , 

i r (1 - f\s{l- sy-^''^F{p, s)G{p, s) =: {F, G) . (A.8) 

Jo J-1 

This follows from the definition of (•, •) by a simple change of variables. Thus instead of the variables 
xi,X2 we have the variables xi = p'^ and X2 = ps. Likewise we write yi = y"^ and 2/2 = yt- Note, 
that in this notation the form of the orthogonal polynomials (jA.ip becomes apparent. 
The operator K,y^,y2 in these variables is given by 

[ICy,tH]{p,s)= [ du'''^{r,h,t2,h)H{E,G) (A.9) 
Ji,J^ 

where 

E = E{p,y;r,h) = (pV + (1 " " + 2py{l - p^)^^! - y^)'/\h)'/^ 

and 

G = G{p,s,y,t : r^MM) = D{C, D{s,t;l,t2);l,h) , 

^ ^ D{p,y;r,ti) D 
E{p,y;r,h) ' E ' 

As before, 

D{a, b; r, t) = ab + {I - a^)^/^{l - b^f'^rt . 
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Recall that 

Kafit) = 1' fiat + sVl^Vr~^)dfi^-'/\s) 
which was used for the Gegenbauer product formula. Now note that 

E J-l 7-1 

The superscript (2) indicates that the operator acts on the second variable of the function. Now 

{F,[ICy,tH]) = 2 C f dmaAr,ti) f dpp^^+\l - p'')" 
Jo J-i Jo 

X ds{l-sy-^'^F{p,s)[Kf^K^^H]{E,s) . 

7-1 E 

Now we proceed using Koornwinder's change of variables: First by going to cartesian coordinates 
we get 



/ / dm^,p{rM)[Kf^K^SH]{E,s) 

Jo J-l E 

/oo poo 
du / dv{l -u"- v%-^-\^^[Kf^K^SH]{E,s) 
-oo Jo E 



where E and D expressed in these new variables are given by 
E 



(1 - - + [(1 - p^)V^il - y^f/^u + py]] '^\d = {1- _ y2)l/2^ + 



Here la,i3 is the normalizing constant. By scaling the variables {u,v) ^ (1 — p^){l — y'^)]^^'^{u,v)we 
get 

' dm^,p{r,h)[Ki^^K^l^H]{E,s) 

J-l E 

OO roD 

2^ ^2 _ ^2^a-/3-l 



/OO f'CO 
du dv{{l-p'){l-y 
-oo Jo 

X v^^[Ki^^K^l^H]{E,s) 



E 



Shifting u ~> u + py and reverting to polar coordinates we obtain 

'1 r'i 



[ f dm^^p{r,ti)[Kf^K^SH]{E,s) 

Jo J-l E 

/•OO /"l 

= - - y^)]- / drr^P^' / da{l - a^f-'l\l - p-^ - y^ - + Ipyraf-^-^ 

Jo J-l 



:(2)^(2) 

I 

where this time 



X [K'^'^K'^>H]{E,s) 

E 



E = r , D = ra . 



CGL September 18, 2009 



32 



Collecting the terms we obtain 

= 2/„,^(l-y2)- [\pp^^+^ 
Jo 

/"OO rl 

X / drr^^+' da{l-ay-^l\l-p'-y^-r^ + 2pyra)%-^-^ 



X / ds{l-sy--'^I^F{p,s)[KPK^^'^H]{r,s) 



-1 



(2) (2) 

The operators and are selfadjoint with respect to the scalar product with the mea- 
sure ds{l — s^)^ . Moreover, they commute which follows from the fact that they have the same 
eigenvectors the ultraspherical polynomials. Hence 

Jo 

poo f'l 

X / drr^'^+i / ^^(i_^2)/3-i/2(i_^2_y2_^2^2pym)°-'^-^ 
Jo J -I 

X I'^dsil - sy~'/^[Ki'^Ki'^F]ip,s)Hir,s) 

which, since the expression is symmetric in p and r, equals {[ICy^tF], H) . ■ 
An similar argument gives the product formula for triangle polynomials first derived by Koorn- 
winder and Schwartz. Recall the scalar product for the orthogonal polynomials on the triangle: 

{f,g):= dxi{l-xi)° {xi-X2)'^x2f{xi,X2)g{xi,X2). (A.IO) 
Jo Jo 

The polynomials, orthogonal in this inner product, are denoted by 



The product formula is 



Krixi,X2) = RTk'\'ixi - l)x\K'''{2^ - 1). 

X\ 



Kf^xlxl)R:i'^iylyl)= [ Rlf^{E\E^H')d^i'^^^''< (A.ll) 

where 

dp"'l^'"'(ri,r2,rs,ri,xpi,ip2,ip3) = dmi3^^{ri,Tp3)dmi3^^{r3,-ip4,)diy^'"'~'^^'^{r2)drna,i3+'y+i{ri,Tpi) , 
with dma<p given by equation ()1.22p . 

^j,/3,7-i/2(^2) = c/3,^(l - r2frl-^'^dr2. (A.12) 

and where 

H = H{xi,X2,X3,Xi;ri,r2,r - 2,r3,ipi,7p2,ip3) 

= E ([l-r2)C^ + r2]'/^E{^,y^;rs,iJ2);rA,i^3) ■ 
V xi yi J 
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The definitions for E, D and C are as in the biangle formula. 

It is convenient to rewrite the inner product in terms of the functions 

F{xi,X2) := f{xl,xlxl) , G{xi,X2) := g{xl,xlxl) , 

so the inner product is given by 



if, 9) 



-. [\xixf^+''^+'\l - xlr f dx2xl^^\l-xlfF{xi,X2)G{xuX2)=: {F,G) . 
Jo Jo 



This fohows from the definition of (•,•) by a simple change of variables. We also let (y 1,2/2) ^ 
(2/1 ) 2/12/2)- "^^^ triangle product formula of Koornwinder and Schwartz is now given in terms of the 
following operator: 



[Gyi.y2G]{xi,X2) = [ [ (imc,,^+^+i/2(n,il) / ^^"^{f 
Jo J-1 Jo 

/ / dmi3^^{rs,t2) / / dmi3^^{r4,t3)G{E,H) 
Jo J-i Jo J-1 



where 

E = E{xi,yi;ri,ti) = {xlyf + (1 - x^il - yfyf + 2xm{l - xi)^/\l - ylf'^nhf'^ 

^ ^ D{xi,yi;ri,ti) 
E{xi,yr,ri,ti)' 

where generally 

D{a, b; r, t) = ab + {1 - a^^/^il - b^f/'^rt . 
Gasper's operator can be rewritten as 

K,Sm = [ j\ fiiy't' + (1 - t'y/'rj + 2yi(l - t')'/\l - yyr2h)r'dm^Ar2M), 
Thus 

[<?o^'^^[Jri')cHHV2,o^K£^,^2) = j\mp,,{r2M f_^dmp,,{nM)G{E,H) . 

As above the superscript (2) indicates that the operator acts on the second variable of the function. 
In the formuls below the constant CQ,,^,-y denotes the products of the various constants normalizing 
the measures that we use. 

Thus the inner product can be written as, 

{FAQyuvM) = ^ fil-xlrxf^"--^^^^ f dvP^^-^/\r2) f f dm.^^+^+i/^Cn, ti) 
JO Jo Jo J-1 

Now we proceed using Koornwinders change of variables and following the discussion of the 
biangle formula we see that E and D become E = r , D = ra. Collecting terms and making the 
final change of varialbes u = [{1 — r2)cr^ + r2]^/^ we obtain 
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Jo Jo J-1 

f du{l - u^iu' - ^2)^-1/2^(1 -xl-yl-r^ + 2x1^1^)^^-^-' 

J 

The proof of Gasper's theorem shows that the operators Ky^^^'^'^ and K^q'^'^^ are selfadjoint with 
respect to the scalar product with respect to the measure dx2{i — X2)^ x'p~^^ . Moreover, they com- 
mute which fohows from the fact that they have the same eigenvectors, i.e., the Jacobi Polynomials. 
Hence the self adjointness follows as in the biangle formula. 

The polynomial preservation also follows from an argument similar to the biangle formula. Here 
we get even powers of E and H and use the fact that the integrals over ti.-.t^ are symmetric. 
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